Finite-temperature dynamical magnetic susceptibility of quasi-one-dimensional 
frustrated spin-1/2 Heisenberg antiferromagnets 
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We study the dynamical response of frustrated, quasi one 
dimensional spin-1/2 Heisenberg antiferromagnets at finite 
temperatures. We allow for the presence of a Dzyaloshinskii- 
Moriya interaction. We concentrate on a model of weakly 
coupled planes of anisotropic triangular lattices. Combin- 
ing exact results for the dynamical response of one dimen- 
sional Heisenberg chains with a Random Phase Approxima- 
tion (RPA) in the frustrated interchain couplings, we calcu- 
late the dynamical susceptibility in the disordered phase. We 
investigate the instability of the disordered phase to the for- 
mation of collective modes. We find a very weak instability 
to the formation of incommensurate magnetic order and de- 
termine the ordering temperature and wave vector. We also 
determine the effects of uniform magnetic fields on the order- 
ing transition. 



I. INTRODUCTION 

A defining property of quasi one dimensional (ID) 
magnets is the weakness of the interchain coupling, whose 
presence nevertheless usually leads to three-dimensional 
ordering. It is natural to consider the ratio of the tran- 
sition temperature Tc to the bandwidth of the excitation 
spectrum D as a quantitative measure of one dimension- 
ality. In systems where Tc/D <C 1 the ordering transition 
occurs in a state where spins on each chain are already 
highly "collectivized" . This means that a very signifi- 
cant proportion of the spectral weight is concentrated in 
a piece of the spin-spin correlation function, which is es- 
sentially one dimensional [1-4]. This suggests that the 
most natural approach to the problem of phase transi- 
tions in such systems is to treat them as instabilities [5], 
driven by weak interchain interactions, of an ensemble of 
ID chains. This approach utilizes the knowledge of the 
correlation functions of individual ID chains, obtained 
by various non-perturbative methods. In this way one 
will automatically reproduce a distinct feature of quasi- 
ID magnets, namely the presence of a broad incoherent 
continuum in the dynamical structure factor. Hence this 
approach is quite different from the conventional spin 
wave theory, which uses individual spins as elementary 
building blocks. Spin wave theory is known to work well 
for the coherent (single-particle) parts of the spectra, but 
it is notoriously difficult to obtain the incoherent parts 
within this approach. Therefore it works poorly in one 
dimension, especially for spin-1/2 magnets, where the in- 



coherence is very strong. The latter is due to the fact that 
excitations of a spin-1/2 Heisenberg chain carry quantum 
numbers different from those of spin waves. 

Some of us have already applied the approach based 
on weakly coupled chains to describe the ordered phase 
of cubic lattice quasi-lD antiferromagnets [6] (see also 
[30,35]). In the present work we apply an analogous 
method to the disordered phase and the ordering tran- 
sition in quasi-lD, frustrated spin-1/2 antiferromagnets 
with Dzyaloshinskii-Moriya (DM) interaction. The ef- 
fects of frustration on quantum magnets are very inter- 
esting and continue to attract much experimental [4,7-11] 
and theoretical [12-18] attention. 

We focus on calculating the dynamical magnetic sus- 
ceptibility xi'-^jk), which is related to the dynamical 
structure factor S{uj, k) relevant for inelastic neutron 
scattering experiments by 
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1 - exp(-a;/T) 



Imx"''(c^,fc). 



(1) 



Note that throughout this paper we will mainly present 
plots of —lTax{ijJ, k) rather than of the structure factor 
itself; the trivial prefactor in (1) can be easily restored if 
desired. 



II. THE MODEL 

Our work is inspired by the recent experiments on the 
frustrated spin-1/2 magnet CS2CUCI4 [4,7-9]. It was sug- 
gested in [4,9] that CS2CUCI4 is described by a Heisenberg 
model on weakly coupled planes of anisotropic triangular 
lattices. Within the planes the exchange couplings are as 
indicated in Fig.l. 



FIG. 1. Exchange paths within the planes: solid lines de- 
note the strong exchange , dotted lines the weaker, frus- 
trated exchange J_l ■ 
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Neighbouring planes arc coupled by a weak antiferro- 
magnetic exchange . In CS2CUCI4 neighbouring planes 
are slightly shifted with respect to one another, but in 
order to keep things simple we ignore this shift. The 
effective magnetic Hamitonian is thus given by 



^ - 2^ ^plane + ^: 



interplane 



(k) ^ 

plane 



o/(fc,fe+l) 

interplane 



Jz ^ ^ Si,j,k-Si,j,k+1 • 



(2) 



In the present work we analyze the model (2) in a 
quasi-lD framework. We choose to view the Hamilto- 
nian (2) as spin- 1/2 chains with exchange J||, coupled 
by a weaker, frustrating, nearest-neighbour exchange J± 
within a plane. Finally there is a very weak antiferro- 
magnetic coupling between neighbouring planes. 

Experimental estimates for the exchange couplings in 
CS2CUCI4 arc J|| = 0.37 mcV, J^/J\i ~ 0.33(1), and 
Jz/J\\ ~ 0.05 [4,19]. Although the interchain coupling 
J±/J\\ is considerable, the smallness of the ratio of tran- 
sition temperature to bandwidth Tc/ttJ\\ « 0.05 indicates 
that the type of approach we are advocating might be ap- 
plicable to CS2CUCI4. 

In addition to the exchange interactions present in (2) 
we allow for the presence of a DM interaction. Our mo- 
tivation is once again the situation in CS2CUCI4, where a 
DM interaction appears to be present [4,7-9,19] although 
it is not straightforward to estimate its magnitude. This 
is because the superexchange between two Cu spin-1/2 
occurs through two Cl~ ligands and not a single one. 
One possible DM interaction that respects the crystal 
symmetry of CS2CUCI4 is 



(3) 



Note that this form of the DM interaction only couples 
spins on a same chain. We mainly concentrate on a DM 
interaction of the type (3) but discuss how to treat more 
general forms in section IX. 

The outline of this paper is as follows. In section III we 
review results for the dynamical susceptibility of a single 
spin 1/2 Heisenberg chain (with exchange anisotropy). In 
section IV we study the effects of a frustrated interchain 
coupling between spin-1/2 chains. We determine the fi- 
nite temperature dynamical structure factor by combin- 
ing the results for single chains with a random phase ap- 
proximation in the interchain couplings. We show that 
an instability towards an incommensurate, ordered state 
develops at a critical temperature, which we determine as 



well. In section V we study the effects of a uniform mag- 
netic field. In sections VI, VII and VIII we carry out the 
analogous programme for spin-1/2 chains with DM inter- 
action. In section IX wc show how to treat more general 
DM interactions within our calculational scheme. Section 
X contains a summary of our results and our conclusions. 



III. DYNAMICAL SUSCEPTIBILITY OF A 
SINGLE SPIN-1/2 XXZ CHAIN 



Let us first consider a single Heisenberg chain with 
"Xy-like" exchange anisotropy < A < 1 

Wxxz = J|| E + ^J^jVl + ^^1^+1- (4) 



We will be particularly interested in the regime A ss 1. 
The spectrum of low-lying excitations consists of a two- 
spinon continuum [20] and is shown in Fig. 2 for the 
isotropic case A = 1. There are gaplcss modes at small 
momentum and at momentum tt, where most of the spec- 
tral weight is concentrated. Spinons carry quantum num- 
ber 5 = ±1/2 and no one-spinon excitations are present. 




FIG. 2. Two-spinon dispersion for the isotropic Heisen- 
berg chain. Also indicated are the sections for which we plot 
the dynamical susceptibility. 

The large distance behaviour of the finite temperature 
dynamical susceptibility can be determined by combining 
results obtained from the Bethc Ansatz solution [21,22] 
of (4) with field theory techniques [23-27] 



x'^'^(T,x) = (-ir^,(77) 



1/2 



sinh (^(a; — iur)) sinh (^(x -|- mr)) 
Here r] is related to the exchange anisotropy A by 



(5) 
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and 



7] = 1 — arccos( A) /n 



u = Jiiao sin(7rry)/(2 — 2t]) 



(6) 



(7) 



is the spin velocity. 

The nonuniversal amphtude Ax{r]) has been deter- 
mined by Lukyanov and Zamolodchikov [28] 



M'n) 



exp 
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°°dt 



sinh(?7i) 



g t Lsinh(t) cosh((l — 77)t) 



77 e 



(8) 



Finally. A is a nonuniversal scale, which recently has been 
calculated in the isotropic case [29] 



A/J|| « 24.27 . 



(9) 



The third factor on the RHS of (5) stems from a 
renormalization-group (RG) improvement in the leading 
irrelevant perturbation to the conformally invariant scal- 
ing limit of (4) [26]. There arc further logarithmic cor- 
rections to (5), some of which have been calculated in the 
isotropic case [27] and lead to (small) corrections to some 
of the formulas wc give below. Wc note that the RG im- 
provement in (5) was done using the inverse temperature 
as the RG scale. It would be interesting to investigate the 
effects of working with two scales i.e. the inverse temper- 
ature and the Euclidean distance. The RG improvements 
are important only in the vicinity of the isotropic point 
A = 1. 

The frequency and momentum dependent dynamical 
susceptibility is obtained by Fourier transformation and 

analytic continuation of the time-ordered imaginary time 
correlation function (5) (see [24,29,31,32]) : 



47rT 



r(i-|-z. 
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AttT 



4ttT 
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Here $(T) is given by 



sin(7rf)r2(l-f) 



1 - 



^\ 4'7-4 

T 



27rT 



(10) 



2-r, 



(11) 



The result (10) is valid for momentum transfers close to 
the antiferromagnetic wave number q ^ ir^ low energies 
a;/J|| <C 1 and low temperatures T / J\\ <C 1. 

The result for the isotropic point A = 1 is obtained 
by taking the limit ^ 1. In this limit diverges and 

y 1 — (A/T)^**"* goes to zero according to 



Ax{ri) 



2(27r)3/2^/T^ ' 



1-(A/T) 



4f;-4 



(12) 

2^r^Vln(A/r) , (13) 



so that their product tends to ^ln(A/r)/(27r)-''/2. 

In Fig. 3 we show the imaginary part of the trans- 
verse susceptibility in the isotropic case for the "constant- 
energy scans" indicated in Fig. 2. The two-spinon contin- 
uum is clearly visible. The maxima at fixed frequency 
occur close to the boundaries of the two-spinon contin- 
uum. 
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FIG. 3. Imaginary part of the dynamical susceptibility 
— Imx'"'(t<J, fe) in units of l/Jy for the isotropic Heisenberg 
chain as a function of q for four different values of w/Jy at a 
temperature of ksT / J\\ = 0.01. 

The longitudinal dynamical susceptibility is given by 



r 1 



(' 



J 

4jj 

' iivT 



4vT 



r M _ 1 — ^' ui+uk \ ' 
\^ 4ri * 47rT J 

where $'(T) is given by [33] 

sin(7r/2r/)r2(l- 1/277) 



(14) 



$'(T) = -AM- 



2nT 



2'?-2n I 



1 + 



A 



2r)-2 



AM = -o 



r(2(i^) 



exp 



L2^r(^)J 

(U / sinh((277 - l)t) 
t \sinh(r7t) cosh((l — 77)^) 



277-1 



-2t 



(15) 
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For A < 1 the longitudinal structure factor at low en- 
ergies is weaker than the transverse one. In Fig. 4 we 
compare the transverse and longitudinal susceptibilities 
for an XXZ chain with A = 0.9. 







m=0.05 


1 


h 





17ii;i6 1571/16 







(0=0.10 







FIG. 4. Minus the imaginary part of the transverse (full 
line) versus longitudinal dynamical susceptibility (dashed 
line) in the anisotropic Heisenberg XXZ chain with A — 0.9. 



IV. DYNAMICAL SUSCEPTIBILITY FOR 
COUPLED HEISENBERG CHAINS 

Let us now take into account the interchain couplings 
J± and Jz- We do this by a Random Phase Approxima- 
tion (RPA) [34,35], which in the present context can be 
understood as the leading term of an expansion in the 
inverse coordination number of the lattice. The dynami- 
cal susceptibility of (three-dimensionally) coupled chains 
are of the form 



l-2J{k,ky,k,)x°"'iio,k) 



(16) 



where a — x,y, z and x°'°'{^^ k) is the dynamical suscep- 
tibility of a single chain. The Fourier transform of the 
interchain spin-spin couplings J{k, ky, k^) is given by 



J[kx,ky,kz) — JzCOs[kz) 

+J± {cos{ky) + cos{kx 



ky)) . (17) 



We recall that we have assumed the coupling along the 
z-direction (between planes) to be antiferromagnetic, i.e. 
unfrustrated. We can recover the results for a weak, frus- 
trated coupling between planes with <^ J± by simply 
setting = in our formulas. Whenever we thus refer 
to = in the following, this should be understood as 
corresponding to a weak, frustrated coupling along the 
z-direction. 

In the experiments on CS2CUCI4 the dynamical struc- 
ture factor was measured for momentum transfers along 



the chain direction. We therefore concentrate on such 
momentum transfers in the framework of our approach. 

The elementary cell and the first Brillouin zone of the 
planar triangular lattice are shown in Fig. 5. If one probes 
the magnetic properties of the crystal along the direction 
of the spin chains, the wave vector transfer follows the 
trajectory ky = kx/2 mod (27r) in the first Brillouin zone 
and has been represented by three arrows on figure 5. 
We see that the wave number transfer along the chain 
direction varies in the interval [0, 47r]. 




FIG. 5. Primitive cell of the triangular lattice and corre- 
sponding first Brillouin zone. Wavenumbers corresponding to 
momentum transfer along the chain direction are indicated 
by arrows. 
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FIG. 6. Imaginary part of the dynamical susceptibility 
— Imx3d(<jJ, k) for momentum transfers along the chain direc- 
tion. The parameters are J±/J\\ ~ 0.1, r/Jy = 0.01 and 
= 0. We see that the frustrated interchain coupling leads 
to an asymmetric line shape. 
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FIG. 7. — Imx3d(<i', fc) for momentum transfers along the 
chain direction. The parameters are J±/J\\ = 0.1, 
T/J|| = 0.02 and J'^/Jn = 0.01. 

In Figs. 6 and 7 we plot (minus) the imaginary part of 
Xsdi'^' ^) momentum transfers along the chain direc- 
tion for four different values of w/J||. Wo have chosen 




k — TT 

FIG. 8. Dynamical structure factor S^^{uj,k) for momen- 
tum transfers along the chain direction. The parameters are 
Jx/J|| = 0.2, T/J|| = 0.05 and JVJy = 0.02. 



parameters such that J^/ J,, = 0.1, T/J|| = 0.01, = ^_ instability condition and critical temperature 



and J_l/J|| = 0.1, T/Jy = 0.02, = 0.01 respectively. 
In both cases the temperature is chosen to be above the 
transition to an ordered state (see below) as it must for 
our approach to apply. 

Figs. 6 and 7 show very clearly that the main ef- 
fect of the frustrated interchain coupling is to make the 
line shape asymmetric. The susceptibility of uncoupled 
chains for momentum transfer along the chain direction 
is symmetric around the antiferromagnetic wave num- 
ber fc = TT as can be seen for example in Fig. 3. This 
symmetry is lost when the interchain coupling is taken 
into account. The asymmetry increases with increasing 
Jj_. The curves for Jj, = and Jz = 0.01 arc qualita- 
tively similar. Another important feature of the imagi- 
nary part of the dynamical susceptibility is that it looks 
incommensurate - the maximum occurs at an incommen- 
surate value of the momentum transfer along the chain 
direction. This is most easily seen in a contour plot such 
as Fig.8, where the dynamical structure factor S'^^(cl), k) 
is shown as a function of energy and momentum trans- 
fer along the chain direction. The dynamical structure 
factor in the disordered phase of CS2CUCI4 displays both 
incommensurabilities and an asymmetric structure [4]. 



The RPA expression for the dynamical susceptibility 
is valid in the disordered phase, i.e. at temperatures 
above the transition to an ordered state. One can ob- 
tain an estimate for the ordering temperature Tc by con- 
sidering the condition for an instability to develop in 
form of a zero-frequency pole in X3d(<^)fc) [34,35]. This 
pole indicates that the system is unstable with respect 
to developing collective modes. In the case of an anti- 
ferromagnet on a cubic lattice the instability signals the 
spontaneous breakdown of spin rotational symmetry via 
the emergence of Neel order and the collective modes 
are the corresponding Goldstone modes (spinwaves). For 
the case of a quasi- ID antiferromagnet on a cubic lattice, 
the critical temperature obtained in this way is compa- 
rable to the experimentally measured Neel temperature 
for KCuFa [35,2]. The leading corrections in the inverse 
coordination nuinbcr of the lattice can be calculated as 
well [36]. In our case the instability condition reads 



2J(fc,fc^,fc,)x""(0,fc) 



= 1, 



T=Tc 



(18) 



where a = x,z. We find that for a = 2; the condition 
(18) leads to a weaker instability, i.e. occurring at a 
lower temperature, except in the case of zero exchange 
anisotropy. We therefore concentrate on a = a; from 
now on. The instability will develop at the maximum 
of J{k,ky,kz)')^^{Q,k). The maximum of the suscepti- 
bility of a single chain x™(0, fc) occurs at = tt, but 
because J{k, fc/2, 0) = vanishes at this point the maxi- 
mum of J')^^ will be shifted away from the antiferromag- 
netic wave number. 
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Extremizing J(fc, fcj,, fc^) x™(0, fc) with respect to ky 
and kz yields ky = k/2 mod (27r) and k^ = n. Setting 
ky = k/2, kz = IT we find that the maximum of the re- 
sulting expression occurs at a value k = ir + ko such that 



7rsinh(Mfco/2Tc 



cos(fco/2) 



cosh(uA;o/2Tc 



cos(^) 



^ , ri .ukn 



u Jz/J± + 2 sm{ko/2) 
0, (19) 



where \l/(x) is the digamma function. Equations (19) and 
(18) with k = n + ko, ky = {tt + ko)/2, kz = tt constitute 
two equations for the two unknowns ko and Tc- Let us 
emphasize that a solution with fco 7^ corresponds to an 
instability at an incommensurate wave number tt + fco- 
This suggests that the resulting order occurs along the 
chain direction with characteristic wave number n + ko- 
In general (18) and (19) can only be solved numerically, 
which is easily done. An exception is the special case 
Jz = and A = 1 (isotropic point), for which we can 
derive analytic expressions. 



1. Isotropic point (rj = 1) 

In the absence of an exchange anisotropy equation (19) 
simplifies further 



uko 



TT tanh(uA;o/2Tc) 



2Im* ( i 



. ukf) 



Solving for xq = fco/(47rTc), we find 
hu ko 



Xo = 



ao 47rTc 



±0.31 , 



= 0. 
(20) 

(21) 



where wc have restored units (fco is measured in units of 
the lattice spacing ao and T in units of ks)- We note 
that fco 7^ 0, which means that the instability develops 

at an incommensurate wave number close to tt. 
The remaining equation (18) then becomes 

2J(fco, 0, 0)x(0, fco) = -2( J. + 2Ji sin(fco/2))x(0, fco) 



1 



(27r)3/2 



T{l+ixo) 



J. 



« 0.25 (^^ + 8 xo ^) Vln(A/Te) = 1 . (22) 

Here we have used the result for the spin velocity of the 
isotropic Heisenberg spin 1/2 chain u = 7rJ||ao/2. 

Let us first ignore the effects of the coupling between 
planes. This case would correspond to a situation where 
the coupling in z direction is also frustrated, but much 
smaller than J±. The critical temperature and ordering 
wave number are then given by 



Te« Aexp(-2.60(J||/Ji)2) , 
|fco| w 2.48^ exp (-2.60 {J\\/J± 

J\\ \ 



(23) 
(24) 



These clearly show that the instability due to the frus- 
trated interchain coupling is extremely weak. Firstly, the 
critical temperature is orders of magnitude lower than its 
counterpart for the square lattice. This is an indication 
how weak the tendency to order is on a frustrated lattice. 
Even for a large interchain coupling J\\/J± = 0.4, the 
critical temperature is only Tc ~ 9.10~® K. The ordering 
wave number is similarly small. Our results should be 
compared to the zero-temperature results for the order- 
ing wave number of [14] and of [15]. The authors of [14] 
used linked-cluster expansion methods to study ground 
state properties and excited states at zero temperature 
of a rather general frustrated model, which contains the 
anisotropic triangular lattice as a special case. They ob- 
tain numerically (in their notations our antiferromagnetic 
wave number corresponds to 7r/2) 



q ~ arccos ( 



1 



V4-6J11/J 



(25) 



which yields a much larger incommcnsuration than our 
result. In [15] a Heisenberg model on an anisotropic 
triangular lattice was studied by enlarging the spin- 
rotational symmetry group from SU(2) to Sp(N) and then 
carrying out a large-N expansion. It was concluded that 
for weak J± < 0.125J|| there is no incommcnsuration at 
all. Given the approximate nature of our results, they 
are certainly compatible with such a scenario. 

In our approach the instability in the SU(2) invari- 
ant case is actually caused by the marginally irrelevant 
current-current interaction that gives rise to logarithmic 
corrections. The weakness of the instability is due to 
the fact that it is caused by precisely these logarithmic 
corrections. 



2. General case 

In the general case 7^ it is not possible to obtain 
simple analytic expressions like (23) and (24) for Tg and 
fco respectively and wc have to resort to a numerical solu- 
tion of (18) and (19). We present the results for various 
values of Jz and J± in Table 1. 
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=0.1 


Jx/J|| =0.2 


Ji/J|| =0.3 





tc = 2.9 10""^ 


= 1.4 10"^^ 


tc = 6.9 10-" 




fco = 7.3 10-112 


ko = 3.6 10-27 


fco = 1.7 10-1° 


0.01 


tc = 0.015 


tc = 0.017 


ic = 0.019 




ko = 0.005 


fco = 0.011 


fco = 0.019 


0.02 


tc = 0.029 


tc = 0.032 


tc = 0.036 




ko = 0.009 


fco = 0.020 


fco = 0.034 


0.04 


ic = 0.056 


tc = 0.060 


tc = 0.066 




ko = 0.017 


ko = 0.036 


fco = 0.060 


0.1 


= 0.13 


= 0.14 


tc = 0.15 




fco = 0.036 


fco = 0.076 


fco = 0.13 



Table 1: Transition temperatures tc = Tc/ J\\ and ordering 
wave numbers fco for isotropic exchange interaction and vari- 
ous values of the frustrated (Jx) and antiferromagnetic (J^) 
interchain couplings. 



susceptibility, ordering temperature and ordering wave 
number. To this end, we must modify the dynamical 
susceptibility of a single chain by taking into account the 
effect of the field. This can be done by utilizing exact 
results obtained by the Bethe Ansatz [22] (see the ap- 
pendix of [37] for a brief summary). For simplicity we 
consider the case of a single, isotropic Heiscnberg XXX 
chain (A = 1). In zero field we have full spin rotational 
symmetry and transverse and longitudinal susceptibilites 
coincide. There are gapless, transverse and longitudi- 
nal excitations at fc = tt and the spectral weight of the 
dynamical structure factor is concentrated in the region 
around fc = tt. Application of a magnetic field breaks the 
symmetry between transverse and longitudinal suscepti- 
bilities. 



We see that that increasing leads to a significant 
increase of the transition temperature and the ordering 
wave number, whereas increasing the frustrated coupling 
J_L mainly increases the value of fco. The interesting 
observation is that a purely antiferromagnetic coupling 
between the chains leads not only to a higher transi- 
tion temperature (as expected), but also to a signifi- 
cantly larger value of the ordering wave number. If we 
take values similar to those for CS2CUCI4 {J±/J « 0.33, 
Jz/J\\ ~ 0.05) we obtain a transition temperature of 
tc = 0.084 and an incommensurablity of fco = 0.083. If 
we take subleading corrections in the RG improvement 
into account [27] these values increase to tc = 0.108 and 
fco = 0.104 respectively. 

These values (which should be taken with a degree of 
caution as J± is not small) are roughly comparable to 
what is observed for CS2CUCI4 [4], where tc » 0.145 and 
fco « 0.186. As is shown in section VII, the presence of a 
DM interaction can further increase tc and fco- 

We may also study the effects of an exchange 
anisotropy A < 1 along the z-direction is spin space (the 
chain Hamiltonian then is given by(4)). The analysis is 
completely analogous the isotropic case. We find that 
the anisotropy enhances both the transition temperature 
and the ordering wave number. For example, for A = 
0.01J|| and J± = 0.1J|| we obtain fco = 0.007, 
O.OlJii and J± = O.lJii we 



0.8. ,L 



tc = 0.023. For A = 
have fco = 0.013, tc 



0.5, J, -- 
= 0.033. 



V. THE EFFECT OF A UNIFORM MAGNETIC 
FIELD 

As we have mentioned in the introduction, the inspira- 
tion for the present work are the observed properties of 
CS2CUCI4. The effects of an applied magnetic field [7,4] 
are particularly intriguing. Let us therefore turn to the 
study of the effects of a uniform field on the dynamical 





" I (b)" " 

FIG. 9. Two-spinon continuum for the isotropic Heiscnberg 
chain in a magnetic field (a) States with AS^ = 1 (b) States 
with AS'' = 0. 

The transverse susceptibility at low energies is dom- 
inated by the contribution from the gapless states at 
fc = TT in Fig. 9(a), whereas the most important contri- 
butions to the longitudinal susceptibility at low energies 
comes from the gapless states at the incommensurate 
wave numbers tt ± S{H) in Fig. 9(b). 

At T = the large-distance asymptotics of the trans- 
verse spin-spin correlation functions in real space is 



{S-{r,y)S^O,0))c,A,{Hy 



(-1) 



v/a.Q 



y + iu{H)T\niH) ' 



(26) 



The magnetic field dependent spin velocity u{H) and ex- 
ponent ri{H) can be calculated to high precision from the 
Bethe Ansatz by solving certain linear integral equations 
[22] . The amplitude (H) is not known analytically but 
has very recently been calculated numerically for several 
values of H using a density matrix renormalization group 
algorithm [38]. 

The asymptotic behaviour (26) holds for distances 
\T+iy/u{H) I large compared to 1 /H. The asymptotic be- 
haviour of the spin-spin correlator at very large distances 
is changed by the application of a magnetic field, but at 
"intermediate" distances the field does not play the same 
role. In order to utilize (26) we therefore need sufficiently 
large magnetic fields. This is the case we are interested 
in. For small fields (26) holds only at extremely large 
distances and it is better to take the field into account in 
RG improved perturbation theory. 
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The dynamical susceptibility can now be determined 
by Fourier transforming (26) and then analytically con- 
tinuing the result to real frequencies. The calculation is 
analogous to the one in zero field [35] and one obtains the 
following result for the low-energy transverse susceptibil- 
ity of a single, isotropic Heisenberg XXX chain (A = 1) 
in a magnetic field 



X''''{uJ,n-k,H) = ^T,H) 



r 



(jn. 

4 



-u{H)k \ 



p /i ri{H) _ ■ ^-uiH)k \ 
\^ 4 ' AttT J 



■+u{H)k \ 
A-kT ) 



l^-^ 4 ' 47rT ) 



(27) 



where |A;| <C 1 and ^{T,H) is given by 

sin(^,?(i?)/2)r2(l-r,(77)/2) 



^{T,H) = -A,{H) 
-u{H) 



u{H) 



1 2-77(ff) 



2'kT 



(28) 



The presently available results for the unknowns in 
(26), (27) are summarized in Table 2. 



M(H) 








u{H)/J\\ao 


0.05 


0.422 


0.121 


0.837 


1.501 


0.10 


0.792 


0.120 


0.782 


1.398 


0.15 


1.109 


0.118 


0.735 


1.259 


0.20 


1.373 


0.117 


0.692 


1.093 


0.25 


1.585 


0.112 


0.653 


0.911 


0.30 


1.748 


0.106 


0.617 


0.721 


0.35 


1.866 


0.095 


0.584 


0.529 


0.40 


1.944 


0.081 


0.554 


0.342 


0.45 


1.987 


0.061 


0.526 


0.165 



Table 2: Magnetization M{H), spin velocity u{H), exponent 
ri{H) and amplitude Ax{H) for different values of the uniform 
magnetic field H. 

As noted above, the result (27) is least reliable at low 
fields as it is obtained by Fourier transforming the large- 
distance asymptotics of the real space correlation func- 
tion. 

Using the results summarized in Table 2 in (27) we ob- 
tain the transverse susceptibility for a single chain. This 
result can now again be combined with an RPA treat- 
ment of the interchain couplings to obtain an expression 
for the susceptibility of the quasi-lD. Following the same 
steps as in the absence of a field we can determine a crit- 
ical temperature and ordering wave number. 
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FIG. 10. Upper two graphs: critical temperature Tc and 
ordering wave number ko extracted from the transverse sus- 
ceptibility as functions of the applied field B = H/J^\ for 
isotropic Heisenberg chains coupled by a frustrated in-plane 
coupling J± — 0.1J|| and antiferromagnetic inter-plane cou- 
pling Jz = 0.01J|| . Lower two graphs: the same for 
J± = 0.2J||. 

In Fig. 10 we plot the ordering temperature and wave 
number as a func;tion of the applied field. We see that at 
first both increase with the field, but Tc{H) eventually 
goes through a maximum and then decreases. This is in 
qualitative agreement with what has been observed for 
CS2CUCI4. If we switch on an exchange anisotropy along 
the z-direction in spin space the qualitative behaviour of 
the Tc and fco as functions of H remain the same. 

So far we have only considered the transverse suscep- 
tibility. In the presence of a field we still have to check 
for possible instabilities in the longitudinal susceptibil- 
ity as the maxima of X^^{^> H) and X™{Q, k, H) occur 
for different values oi k. It is important to note that 
the physical nature of a longitudinal instability is very 
different from that of a transverse one. The latter is as- 
sociated with the spontaneous breakdown of the spin ro- 
tational symmetry around the direction of the magnetic 
field, whereas a longitudinal instability does not break 
any continuous symmetries, but rather corresponds to 
the formation of a spin-density wave in the ground state. 

The longitudinal susceptibility in a field can be calcu- 
lated along the same lines as outlined above for the trans- 
verse ones. However, the longitudinal real-space correla- 
tion function decays faster with distance as the field is 
increased. Therefore it is likely that Fourier transform- 
ing the large-distance asymptotics leads to a less reliable 
result than for the transverse susceptibilities. It would 
be very interesting to investigate this issue by numerical 
methods. Fourier transformation of the large-distance 
asymptotics and analytical continuation yields 



X 



'{u,n±S-k,H) = ^'{T,H) 
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p ( 1 _ ■ u:-u{H)k \ p ( 1 _ ■ oj+u(H)k \ 



1 _ ■ uj-u{H)k 
iri(H) * 47rT 



)^( 



1 1 ■ i^+u(H)k \ ' 

ir]{H) ' 47rT / 



(29) 



where 5 = 2ttM{H) and $'(T, if) is given by 

sin(V27?(i?))r2(l-l/2r,(i?)) 



$'(T,ir) = -A4ir)- 



27rr 



2-i/n{H) 



(30) 



The amplitude (H) is again known numerically [38] , so 
that we can repeat the analysis carried out above for the 
transverse susceptibility. We find that in general an in- 
stability in the longitudinal susceptibility exists at a crit- 
ical temperature comparable to the one for the transverse 
instability. The corresponding characteristic wave num- 
ber fco is located in the vicinity of tt ± ^ = vr ± 2ttM{H), 
where M{H) is the magnetization of uncoupled chains. 
Note that tt ± 5 are the points where incommensurable 
soft modes exist in the AS*^ = sector of the excita- 
tion spectrum of individual chains (see Fig. 9(b)). For 
large values of the applied field H,-k±5 approach and 
27r respectively. We note that in this regime the further 
complication arises that it is no longer appropriate to ne- 
glect the contributions (as has been done in (29)) of the 
smooth, non-oscillating terms in the spin-spin correlation 
functions. 

As a function of the applied field i?, the transition 
temperature of the longitudinal instability ij°"^(ii) first 
increases, goes through a maximum and then decreases 
again. This is similar to what we obtained for the trans- 
verse instability. However, in contrast to the transverse 
instability the maximum of tj,°"^(-ff) occurs at small mag- 
netization. 

Which instability is dominant is rather sensitive to the 
values of J_l, and the applied field H. For = 0.1 J||, 
Jz = 0.01J|| we find that the transverse instability occurs 
at a higher temperature except for small fields, where the 
longitudinal susceptibility appears to dominate. On the 
other hand, as we have pointed out above, for small fields 
our results for the Fourier transforms are least reliable. 

So far we have worked with a spin rotationally sym- 
metric Hamiltonian. Let us now consider the effects of an 
exchange anisotropy. We have to distinguish two cases, 
depending whether or not the applied field is along the 
direction of the exchange anisotropy. In the former case 
the analysis is completely analogous to the spin rotation- 
ally symmetric case and we find a behaviour Tc{H) and 
ko{H) very similar to Fig. 10. The transition temperature 
increases with increasing field, goes through a maximum 
and eventually decreases again. 

Let us now turn to the case where we have an exchange 
anisotropy in z-direction and apply the field along the 
x-direction. The effect of the magnetic field is now to 



generate an excitation gap. This can be seen as follows. 
The chain Hamiltonian is 

Wxxz,H = -^11 S^S^j^i + SjSj_^_i + ASjSj_^_l 



+hY,s^ = no + ni , (31) 

j 

where Tio is the Hamiltonian of the anisotropic spin-1/2 
chain and 



Hi 



(32) 



Let us study the effect of the perturbing operator by 
bosonizing at the critical point defined by Ho and then 
perturbing this fixed point theory by (32). The bosonized 
form of Ho is 



Ho 



\ j dx[{dM 



(33) 



where $ is a canonical bosonic field and 6 is the dual 
field. The uniform part of the x-componcnt of the mag- 
netization is given by the product of two operators 

H cos {y2^ cos {y2^ $ j , (34) 

where 7] has been defined in Eq.(6). These operators 

are formally relevant and not being Lorentz scalars they 
belong to the class of perturbations with non-zero con- 
formal spin. Such operators require a special treatment 
(see [39], chapters 8 and 20 and references therein). In 
particular, to second order, the following two spin-zero 
fields are generated 



cos(-v/87r/77$). 



(35) 



The resulting problem has been considered in detail in 
Chapter 20 of [39]. Because we have < 1, the renor- 
malization group equations flow to a Sine-Gordon model 
in the dual field. For weak fields H, a spectral gap of 
order of 



MocH 



(36) 



is generated. The gap grows with increasing magnetic 
field. Clearly, the growth of the gap will lead to a de- 
crease in the transition temperature Tc - the effect ob- 
served in CS2CUCI4 when the field is applied along the c 
direction. 

There is a second way to see that application of a uni- 
form field at an angle to the exchange anisotropy induces 
a spectral gap [40]. If we consider the Hamitonian 



WzXX,H = J\\ ^ SjSj_^_i 

j 



QZ QZ 1 A QX QX 



(37) 
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and bosonize the isotropic Hciscnbcrg chain in a field first 
and then take the exchange anisotropy into account as a 
perturbation, we obtain a Sine-Gordon model for the dual 
field [40,37]. The cosine term in the Sine-Gordon model 
is relevant and generates a spectral gap. The spectrum 
of the Sine-Gordon model consists of soliton and anti- 
soliton only and the dynamical structure factor (in the 
xy plane) displays an incoherent two-particle continuum. 
This suggests that the elementary excitations are massive 
spinons. 

VI. A SINGLE SPIN-1/2 XXX CHAIN WITH DM 
INTERACTION 

As discussed in [9] , DM interactions may play a role in 
accounting for all the observed properties of CS2CUCI4. 
We therefore consider now a single, isotropic spin- 1/2 
Heiscnberg chain with a DM interaction along the z di- 
rection in spin space 

L 

'U 7' \ ^ QX QX j_ ay CV I QZ 

(38) 

i 

It is well known that (38) with periodic boundary condi- 
tions is equivalent to an XXZ chain with twisted bound- 
ary conditions. Indeed, a local rotation around the z-axis 

S] = 5| , (39) 
with 9 = — arctan(£)/J') maps the Hamiltonian (38) onto 

L 

Here J = J' /cos 9 and A = cos^. For a system with 
open boundary conditions there are no further changes. 
In particular, we can rest assured that bulk correlation 
functions in (38) can be obtained from bulk correlators 
in the anisotropic spin- 1/2 Heisenberg chain. The ex- 
citation spectrum of a spin- 1/2 Heisenberg chain with 
DM interaction can be obtained from that of the corre- 
sponding XXZ chain by taking into account the shift in 
momentum induced by the mapping (39). In Fig. 11 we 
show the qualitative excitation spectrum with quantum 
numbers AiS^ = ±1 around momentum tt. There are 
incommensurate soft modes at tt ± 0. Excitations with 
AS^ = stay commensurate i.e. become soft at momen- 
tum TT. 




k 

FIG. 11. Schematic two-spinon dispersion in the vicinity 

of fc = TT in the sector AS" = ±1 for the isotropic Heisenberg 
chain with DM interaction. 

Using the above mapping we can now express bulk cor- 
relation functions of the spin- 1/2 chain with DM inter- 
action in terms of correlation functions of an XXZ chain 
with exchange J/ cos 9 and anisotropy A = cos 9. For 

example, 

(^-^;+fe)DM = e-'"' (5r5++,)xxz, (41) 

where wc have used that by global spin rotational sym- 
metry (of the bulk) around the x-axis 

i^j^j+k) XXZ = {Sj Sj_^_^)yixZ , (5'JS'|^j,)xxZ = 0. (42) 

This allows us to express the dynamical magnetic suscep- 
tibility of the model (38) in terms of the results for the 
Heisenberg XXZ chain 

XuMi^,k) = x'^iuJ,k + 9), 

X'^Ui^,k) = x'%oj,k). (43) 

In Fig. 12 we plot — Imx^^^f (cj, fc) for a DM angle of 
6* = 0.1 as a function of q for four different values of 
uj/J\\ at a temperature of A;bT/J|| = 0.01. Had we plotted 
instead 

x''''{u;,k) = ^{x+-{u^,k-9)+x-+{co,k + 9)) , (44) 

we would have seen an incommensurate four-peak struc- 
ture. Because of the chirality introduced by the interac- 
tion, xim s-'^d Xdm iio'w differ (u^.S^k) x S{—k) has a 
nonvanishing expectation value). 
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FIG. 12. Imaginary part of the dynamical susceptibility 
— ImxoM A;) (full line) and —ImxEmi'^jk) (dashed-dotted 
line) as functions of q for four different values of a// J at a 
temperature of ksT/ J = 0.01. The DM angle is chosen to be 
6* = 0.1. 
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FIG. 13. — Imxjjd (solid line) and —\v[ix^^{y>,k) 

(dashed-dotted line) for a DM angle 6 = 0.1, Jx/Jy = 0.1, 
T/ J|| = 0.02 and = 0. The frustrated interchain coupling 
breaks the symmetry around k = t:. 



VII. DYNAMICAL SUSCEPTIBILITY FOR 
COUPLED HEISENBERG CHAINS WITH DM 
INTERACTION 

In the presence of a DM interaction along the z- 
direction in spin space, the RPA result for the dynamical 
susceptibility of coupled chains is given by 

+-( M_ x+-i^,k-e) 

X3/(w, k) = xtdi^' ^' ^ ^ ' 

y^^ (iv k) = X''{^,k) 

' > l-2J{k,ky,k,)x''{uJ,k) ' 



(45) 



where J(k,ky,kz) is given by (17) and x'^'^{'^,k) is the 
dynamical susceptibility of a single Heisenberg chain. We 
note that the RPA expressions for X3^(w,fc), a = x,y, 
arc not simply (16) with x""(t^, k) replaced by x""j(c(;. A;). 
Instead, Xsd 
xtd and X3d 



xld and 



XV 

Xsd = 



are obtained from 



XS = I (XsV + X3d+) ' X3^ = ^ ixti - X3d+) • (46) 

In Figs 13 and 14 we plot the imaginary parts of the dy- 
namical susceptibilities Xsd" i^' Xdm('^' ^) as functions 
of the momentum transfer along the chain direction for 
a DM angle of ^ = 0.1, = (the behaviour for 7^ 
is qualitatively the same) and two different values of J±_ . 
The effect of the frustrated interchain coupling is again 
to remove the symmetry around k = it. 
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FIG. 14. — ImXsd (^^7^) (solid line) and —lmx3^{u!,k) 
(dashed-dotted line) for a DM angle oi = 0.1, Jx/Jy = 0.2, 
T/J|| = 0.04 and = 0. 

Let us now determine the transition temperature and 
ordering wave number in the presence of a DM interac- 
tion. We will see that the effect of the DM interaction is 
to significantly increase both quantities. The instability 
conditions read 



J{k,ky,k,)x+-i0,k-9) = 1 , 

j{k,ky,k,)x-+{o,k + e) = i , 



(47) 



where the ID dynamical susceptibilities fulfil x^ (0, k) = 
X ''(0, k). As before the instability develops at the max- 
imum of J{k,ky,kz)x~^~{0,k ± 9). In the absence of 
a DM interaction {9 = 0) this maximum was shifted 
away from the maximum of the susceptibility for a single 
chain at A; = tt, because J{k,k/2,n) vanished precisely 
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at fc = TT. In the presence of a DM interaction the max- 
ima of the single-chain susceptibihty occurs at tt ± and 
the effect of the frustrated interchain coupHng is there- 
fore different. We find that the numerical value of k at 
which the instability develops is largely determined by 
the DM interaction and to a lesser degree by the frus- 
tration. However, the frustrated interchain coupling de- 
stroys the symmetry in momentum space around k = tt 
and determines whether the instability will develop in the 
vicinity of tt -|- 6' or of tt — f?. The fact that the transition 
temperature is increased by the DM interaction can be 
understood by considering the susceptibility of a single 
chain with DM interaction. As we have shown in sec- 
tion VI the corresponding Hamiltonian maps onto an an 
anisotropic Heisenberg XXZ chain. It follows from (10) 
that for A < 1 the transverse susceptibility is enhanced 
as compared to the isotropic case A = 1, which in turn 
leads to a higher Tc. 

Solving (47) and the equation for the maximum of 
J{k, ky,kz) X"' (0, k ± 0) numerically, we obtain the re- 
sults for the transition temperature Tc and ordering wave 
number fco shown in Table 3. 
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J±/J\\ =0.1 


Jj_/J|l =0.2 


Jx/J|| =0.3 


0.05 


tc = 0.008 
fco = 0.053 


tc = 0.017 
fco = 0.061 


tc ^ 0.028 
fco = 0.076 


0.10 


tc = 0.016 
fco = 0.10 


tc = 0.032 

fco = 0.12 


tc = 0.051 

fco = 0.14 


0.15 


tc = 0.023 
fco = 0.16 


tc = 0.046 
fco = 0.18 


tc = 0.073 
fco = 0.21 


0.20 


tc = 0.031 
fco = 0.21 


tc = 0.060 
fco = 0.23 


tc = 0.095 
fco = 0.27 



Table 3: Transition temperatures tc = Tc/ J\\ and ordering 
wave numbers fco for various values of the frustrated {J±) 
interchain coupling and DM angle 6. The interplane coupling 
Jz is taken to be zero. The effect of a small Jz on Tc, fco is 
negligible compared to the effect of the DM interaction. 

We see that a strong DM interaction leads to much 
larger values for Tc and fco than a frustrated interchain 
coupling alone. For the values of couplings observed in 
CS2CUCI4 J± = 0.33 J, Jz ~ 0.05 J, D « 0.05 J we obtain 
tc = 0.11, fco = 0.14, which are close to the measured 
values. However, recent evidence suggests that the DM 
interaction in CS2CUCI4 is not of the kind considered in 
this section [19]. 

VIII. DM INTERACTION AND A MAGNETIC 
FIELD 

Finally, let us investigate the case where we have both 
a DM interaction and a magnetic field. Now it is crucial 
whether or not the field is applied along the direction 
singled out by the DM interaction. If it is, the system 



remains gapless and wc can proceed along the same lines 
as before. If the magnetic field is applied at an angle 
to the DM interaction, we believe that an excitation gap 
is generated [42]. Let us consider the former case. The 
chain Hamiltonian is of the form 

L 

-iv T' \ ^ ox ox , ay ay 1 

j j 

By means of the unitary transformation (39) this maps 
onto an anisotropic Heisenberg XXZ chain in a field 

L 

+HJ2S], (49) 
j 

wore J = J'/cos9 and A = cos9 (9 ~ — arctan(T'/J') 
as before). The model (49) remains gapless and wc can 
determine the finite-temperature dynamical susceptibil- 
ity by the same methods wc used in the absence of a 
DM interaction in section V above. The result for the 
transverse susceptibility is of the form (27), where the 
exponent ri{H, D) and velocity u{H, D) can again be de- 
termined from the Bethe Ansatz, whereas the amplitude 
Ax{H) is known numerically [38]. Taking the interchain 
couplings into account in RPA and then looking for insta- 
bilities as before, we obtain the results shown in Fig. 15. 




0.5 1 1.5 0.5 1 1.5 

B B 

FIG. 15. Critical temperature Tc and ordering wave num- 
ber fc() as functions of the applied field B = H / for isotropic 
Heisenberg chains with a DM interaction, coupled by a frus- 
trated in-plane coupling J±_ = 0.1J|| and no antiferromagnetic 
inter-plane coupling = 0. The DM angle is indicated in the 
figures. 

The behaviour of Tc{H) and ko{H) as functions of the 
applied field is similar to what we found in the absence 
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of a DM interaction. However, like for the zero field case 
a strong DM interaction leads to a significant increase in 
the absolute values of Tc and fco- 



IX. GENERAL DZYALOSHINSKII-MORIYA 
INTERACTION 

In this section we discuss how to treat more general 
types of DM interactions within the framework of our 
coupled chain approach. 



A. Dzyaloshinskii-Morya interaction along the chains 

In a compound such as CS2CUCI4, the three dimen- 
sional elementary cell contains four Cu^+ ions [41]. Two 
of them, say (1) and (2), lie within a planar layer and 
are coupled by a relatively strong, frustrated exchange 
interaction. The other two, say (3) and (4), lie in the 
layer above and are coupled weakly to (1) and (2) [9]. 

In the case where it involves spins along the chain di- 
rection, the most general form of the DM interaction 
which is allowed by symmetries [9] is 



n,jn 



(50) 



Generically, this interaction requires the distinction be- 
tween the sets (1) and (2). (Up to now, we had only 
considered the case Di = D2). 

As a consequence, and from now on, the Heisenberg 
chains forming up the triangular planar lattice will not 
necessarily be considered all equivalent, but belonging 
alternatively to sets (1) or (2). This doubles the primitive 
cell chosen in figure 5 in the y-direction. The magnetic 
Hamiltonian that takes the difference between the two 
types of spins into account is 

+ Jj_ ^ '^iX+i • (-^i+l.m + ^n-l,m+\) 



~1~ ^DK 



(51) 



Due to the doubling in the y-direction, the Fourier trans- 
form of the interchain spin-spin couplings is now 



J{k) = 2J± ( cos{ky/2) + cos{kx - ky/2)) 



(52) 



In order to choose the quantization axes for sets (1) 
and (2) we introduce the unit vectors 



ui =Lii/|Z)i|, U2 = D2/\Di\ 

V = U\X U2/\ui X U2I • 



(53) 
(54) 



The spins along the chains of type (1) and (2) will be 
quantized in local coordinate systems with axes {v, U\ x 
v,u\) and {v,U2 x v^u^) respectively. We denote the 
angle between ui and U2 by C, and define two DM angles 
by 0\ = — arctan(|i5i|/J||), 62 = — arctan(|i52|/J||). 

Specializing to CS2CUCI4, we note that the space sym- 
metry group for this material is Prima. This group con- 
tains symmetry elements which map the spin chains of 
different types onto each other, which leads to the fol- 
lowing restrictions of the DM vectors Di and 1)2 • They 
are perpendicular to the chain direction, perpendicular to 
each other and of equal length. This yields that 61 = 62 
in this case; sec [9] for a detailed symmetry analysis. 

Next we define the usual step operators for j = 1 , 2 : 



(55) 



Notice that because of our choice of the quantization axis, 
the letters x, y, z refer to difi^erent directions for j = 1 and 
j = 2. In Fourier space, the total Hamiltonian is 



COS(fcj; — 6j) 

COS 6i 



sl-Ak)S7-,J-k) 



k;j=l,2 

+J|| COS k,Sl^^{k)S(^^{-k) 

+ cos C E •^(^) ['^fi) (^)^f2) i-k) + {k)Sl2) i-k) 

k 

+ Sin C E Ak) [5fi) {k)Sf,^ i-k) - (^)5f2) {-k) 

k 

+ Y.J{k)Sl,){k)Sl2){-k). 



(56) 



Next we will write the "effective" quadratic spin 
Hamiltonian corresponding to the random phase approx- 
imation. In order to do so, we define 



A, 



2 
1 

2i 



{[xnk, + 0,)]-'-[xf{k,-9,)]-'), 



= ixTik.)]-' . 



(57) 



Here x^'^ and Xj^ denote the time-ordered imaginary-time 
correlation functions for chain (j) in the presence of the 
DM interaction with coupling Dj. They are related to 
the transverse and longitudinal dynamical susceptibilities 
for chain (j) by a Wick rotation and a global minus sign. 

The effective RPA correlation functions between spin 
operators (Sf , S'f, S'f , 5*1 , 5*1, S'|)(fc) and spin operators 
{Sf,Sf,S^,S^,S^,S^y{-k) are then given by the in- 
verse of the bilinear form 
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Si 


Ai 





J 








-Ai 


Si 








J COS Q 


J sin 








f7i 





— JsmC, 


J cos 


J 










A2 








JcosC 


—JsinC, 


-A2 


S2 








Jsin^ 


JcosC 








O2 



(58) 

The emergence of a pole in the dynamical susceptibility 
corresponds to the vanishing of the determinant of this 
matrix. This leads to the following necessary condition 

= J6 - [(1 + cos^ C)SiS2 + sin2 C(Si02 + S2O1) 

+ ^11^2 cos^ C - 2Ai A2 cos C] 

+ [(l + COs2C)Oil72SiI]2 

+ sin2 C (r!iSi(S2 + A2) + ^^2^2(^2 + A^)) 

+ cos^ C(Si + A2)(S2 + A^) - 2f]il72 Ai A2 COS C] 



Oifi2(S? + A?)(I]2 + A2). 



(59) 



In the particular case when ui = U2, but 0i is not neces- 
sarily equals to 62, the condition simplifies to 

= [1 - J\T{k)xTm^ - J\T{k + 0i)xr{k + 02)] 

x[l-j\r{k-ei)xr{k-02)], (60) 

where k stands for kx- If in addition to ui = U2 we also 
have 01 = 62 then we recover, as we must, the instability 
conditions discussed before. 

The dynamical susceptibilities can be calculated, by 
inverting (58), then continuing analytically on the fre- 
quencies. 



B. Dzyaloshinskii-Moriya interaction along the 
interchain bonds 

Up to this point we have only considered the case when 
the DM interaction involves spin along the chains. It was 
recently suggested to us by R. Coldea and A. Tennant [19] 
that the DM interaction in CS2CUCI4 may also involve 
pairs of spin along the interchain couplings. 

The direction of the DM vector in this case is again 
constrained by crystal symmetries [19]. The DM vector 
is perpendicular to the triangular planes and is staggered 
along the chain direction as shown in figure (16). 




DM vector pointing downward (•) DM vector pointing upward 

FIG. 16. Schematic representation of the orientation of the 
allowed directions of the DM vector within one layer for the 
case where the interaction involves spins along the interchain 
bonds. 

For CS2CUCI4 the DM vector appears to be staggered 
between triangular layers as well [19]. In what follows we 
will first consider the simpler case in which there is no 
staggering between layers in order to keep our formulas 
simple. 

We again need to distinguish between spin chains of 
type (1) and (2). The calculations are similar to the ones 
of subsection IX A. Since the longitudinal spin operators 
Sq-. do not couple to the transverse ones, the transverse 
and longitudinal susceptibilities can be calculated sepa- 
rately. In Fourier space, the Hamiltonian takes the form 

'^DM = 'Ho + Wint , 

^o = J|| E JZ^oskx S^^){k)S^^^i-k) 

k;j=l,2 " 
Hint=535(^)^ A{k) SH-k) 



+ 5^2J(fc)5fi)(fc)5f2)(-fc) 

fe 

+ E E COS ^0) (^)^O) (-^) ' 

where D denotes the DM coupling, 
Sikf = {S+^{k),S+^{k)) , 
J{k) = J^(cos(fcy/2) + cos{kx — ky/2)), 
K{k) = D{ sm{ky/2) + sin(fc^ - ky/2) + ) , 

and A{k) is the 2x2 matrix 

Jz cos kz J + K 
J + K Jz cos kz 

It will be convenient to define 



(61) 



(62) 



(63) 



Le{k) = J{k) + K{k) 

= J^[cos(^ 
cos0^ ^2 



■ 6) + cos(fc^ 



0)], (64) 



where 6 = — arctan(D/J_L). 
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Wc denote by x"* (k) = x. ^(^) the time-ordered 
imaginary-time transverse correlation function between 
spins S'^j){k) and spin S^j^{—k) in the absence of the in- 
terchain couphng and DM interaction. This is of course 
simply the correlation function of a single one dimen- 
sional chain. 

The RPA expressions for the transverse correlation 
functions between spin operators Sa{k) and iS^(— fc) 
{a, (3 = 1,2) are given by the matrix elements of the 
inverse of the following 2x2 matrix 



[X"*" ] ^ + Jz cos kz Le 

Le [X^"]~ +JzCoskz 



(65) 



The time ordered imaginary time two point correlation 
function of spins is obtained by adding the contributions 
from the various sublattice correlators, i.e. by taking e.g. 



(66) 



After analytic continuation we obtain the following 
RPA formula for the transverse dynamical susceptibili- 
ties 

+-/ A X+"(w,fcx) 

Xj; [UJ, k) = =; , 

1 - [Leik) + J^cosA;^] x+"(w,fca;) 
X3j'{o^,k)=xtdi^^-k)- (67) 

Here x^ {'-^i k^) is the dynamical susceptibility of a sin- 
gle one dimensional chain. We note that we could have 
arrived at this result by first removing the DM interac- 
tion along the interchain bonds by moans of a unitary- 
transformation in a way analogous to (39). This induces 
an effective DM interaction along the chains and the re- 
sulting Hamiltonian can again be analyzed within RPA. 
After undoing the unitary transformation one recovers 
(67). 

The structure of the dynamical susceptibilities (67) is 
quite similar to the one we obtained in the case where the 
DM interaction is along the chain direction (45). There 
again is an effective global shift of the momentum by a 
"DM angle" 9. However, now the angle depends on the 
ratio of the DM coupling D to the interchain coupling 
J± rather than the strong exchange J|| , which makes it 
much larger. This leads to an enhancement of the effect 
of the DM interaction. 

The condition for the emergence of a transverse insta- 
bility towards the formation of a coherent mode reads 

1 = max |[J:e(fc) + Jz coskz] X"^~(w, k^)^ . (68) 

One finds that the instability again occurs along the chain 
direction, i.e. ky = kx- 



C. Taking into account the staggered layers 

It was suggested in [19] that in CS2CUCI4 the direction 
of the DM vector alternates between neighbouring layers. 
This feature can be accomodated within our calculation 
as follows. 

Let us first suppose that the intcrlaycr coupling con- 
stant is negligible. Then the global transverse dy- 
namical susceptibilities are given by the sums over the 
contributions from the two typos of layers. In this case 
the transverse susceptibility reads 

X+-{uJ,kx) (1- J(fc)x+-(w,fc.)) 
(l-Le{k)x+-{oJ,kx)) (l-L_e(fc)x+-('^,fcx)) ' 
X3/('^.^)=X3cr(^.^)- (69) 

As expected the transverse dynamical susceptibility has 
lost its chiral nature, i.e. we now have Xsd^ — xtd ■ 
This is in marked contrast to the case where the DM 
interaction involves spins along the chains, for which we 
always had X3/ 7^ xtd- 

If Jz is not zero the calculations are slightly more com- 
plicated. One has now to distinguish spins in neighbour- 
ing layers and the elementary cell is doubled in the z- 
diroction. We will denote the spin operators correspond- 
ing to the four sites per unit cell by S'g ^.^ (j, k = 1,2). Let 

us also introduce the function I{k) = J^ cos{kz/2). Then 

the effective RPA transverse correlation hmctions be- 
tween spin operators {S^i ly 1)^ 2)^ '^(2 2) ')(k) and 

{S^^ 1); 'S'(2 1)^^(1 2)' '^(2 2))'^(~k) given by the inverse 
of the following bilinear form 



ix' 



Le 
I 




Le I 

[x+-]-^ 

[X+-] 

/ L-, 



L. 

[x+- 



(70) 



The time ordered imaginary time two point function 
of spin operators is again obtained by summing over the 
sublattice contributions. After analytic continuation, we 
obtain the following RPA expression for the transverse 
dynamical susceptibility 



xtd (^' k) = 

X+-{uJ,k,) (l + Niik) x+-{uJ,k,)) 

[l-2J{k)x+-{io,kx)+N2{k) [x+-{^,kx)f) ' 
where 



(71) 



Ni{k) 

N2{k) 



I{k) - J{k) , 
J2(fc) - K^{k) 



I\k). 



(72) 
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We again have that 

Xs^{u;,k) = xtii^,k). (73) 

Prom (71) we obtain a modified set of instabiUty condi- 
tions 

(^J{k) ± ^K^{k)+P{k)^ x+" (0, fcx) = 1 . (74) 

Extremizing with respect to kz and ky yields the condi- 
tions ky = kx and kz = 0. 

If we specify the exchange couplings according to the 
values suggested for CS2CUCI4 (Jz — 0.05J||, J± ~ 
0.33J||, D ~ 0.05J||) we obtain a critical temperature of 
Tc = 0.727 K and an ordering wave number of ko = 0.154. 
These arc close to the experimental values Tc = 0.62 K 
and ko = 0.186. 



6 


Ji/J|| =0.1 


Ji/,/|l =0.2 


J±/J 


1=0.3 


0.05 


tr = 0.070 


t,, = 0.078 


tc = 


0.092 




A'o = 0.021 


Ay = O.U4() 


/U) = 


0.082 


0.10 


tc = 0.074 


tc = 0.091 


tc = 


0.120 




ko = 0.022 


ko = 0.053 


ko = 


0.103 


0.15 


tc = 0.079 


tc = 0.109 


tc = 


0.153 




ko = 0.023 


ko = 0.062 


ko = 


0.128 


0.20 


tc = 0.087 


tc = 0.130 


tc = 


0.189 




ko = 0.025 


ko = 0.073 


ko = 


0.153 



Table 4: Transition temperatures tc = Tc/J\\ and order- 
ing wave numbers ko for various values of the frustrated 
(Jx) interchain coupling and interchain DM angle 6 = 
arctan(_D/Jx). The interplane coupling Jz is taken to be 
0.05 J|| . 

X. SUMMARY AND CONCLUSIONS 

In this work we have studied the dynamical response of 
frustrated quasi-lD spin-1/2 Heisenberg magnets in the 
disordered phase. The starting point of our approach arc 
exact results for the finite-temperature dynamical suscep- 
tibility of an ensemble of uncoupled chains. Taking the 
couplings between chains into account within the frame- 
work of a random phase approximation, we obtained an 
analytic expression for the dynamical structure factor of 
the quasi- ID system we are interested in. In the dis- 
ordered phase the main effects of the frustration are to 
generate an asymmetry of the line shape and a shift of the 
apparent dispersion to an incommensurate wave vector. 

By analyzing the instability of the disordered phase 
with respect to the formation of collective modes wc stud- 
ied the transition to the low-temperature ordered phase. 
In particular we determined the ordering temperature 
and ordering wave vector within the framework of our 
approach. We found that there is a very weak instability 
towards an incommensurately ordered state. 



We then considered the effects of an applied magnetic 
field. We found that for isotropic Heisenberg magnets 
application of a magnetic field leads first to an increase 
in the transition temperature Tc{H) and for very large 
field to an eventual decrease. In the presence of an ex- 
change anisotropy we found two distinct behaviours: if 
the field is applied along the direction of the anisotropy 
the situation is very similar to the isotropic case. Ap- 
plying the field at an angle to the anisotropy generates 
a gap in the individual chains, which leads to a decrease 
of Tc{H) with H in the quasi- ID system. This type of 
magnetic phase diagram is qualitatively similar to what 
has been observed in CS2CUCI4. 

In the second part of this work we took into account 
the effects of various types of Dzyaloshinskii-Moriya in- 
teractions, which break spin rotational invariancc. We 
first considered DM interactions involving spins along the 
chain directions and then DM interactions involving spins 
along interchain bonds. In the disordered phase the main 
effect of DM interactions is to generate stronger incom- 
mensurations. DM interactions also leads to a significant 
enhancement in the transition temperature and ordering 
wave number. 

In general our results are qualitatively similar to what 
is observed experimentally in CS2CUCI4. In presence 
of a DM interaction of the kind proposed in [19] for 
CS2CUCI4 we obtain a transition temperature and or- 
dering wave number, which arc close to the experimental 
values. However, given that the interchain coupling is 
not small it is unclear how reliable the RPA approach 
is. In light of this fact it would be very interesting to 
determine the leading corrections to RPA, which can be 
done for example along the lines of [36]. It also would 
be interesting to extend the coupled chains approach to 
the ordered phase. This is not straightforward due to the 
presence of "twist" -operators [17]. 
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